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CRUSTAL  AND  UPPER  MANTLE  STRUCTURE  OF  QINGHAI-XIZANG  PLATEAU 
AND  ITS  SURROUNDING  AREA 


Ta-liang  Teng 

(Department  of  Geological  Sciences, 
University  of  Southern  California, 
Los  Angeles,  California,  USA  90007) 


ABSTRACT 

Multiple-filter  technique  and  generalized  linear  inversion  are  applied 
to  group  velocity  data  derived  from  properly  rotated  three-component  digital 
surface  waves  from  Seismological  Research  Observatory  stations.  Repeated 
group  velocity  measurements  over  the  same  paths  were  made  possible  by  using 
large  aftershocks  of  several  great  (M  >  7)  earthquakes  that  recently  occurred 
in  China.  Specific  attention  is  paid  to  the  crustal  and  upper  mantle  struc¬ 
ture  of  the  Qinghai-Xizang  plateau  and  its  comparison  with  those  of  neigh¬ 
boring  tectonic  provinces.  The  results  demonstrate  the  strong  lateral 
heterogeneity  in  the  crust  and  upper  mantle  of  the  Chinese  mainland.  A 
particularly  thick  crust  and  an  upper  mantle  of  relatively  low  velocity  are 
found  underlying  the  Qinghai-Xizang  plateau  that  corroborates  the  continental 
collision  between  India  and  Eurasia. 


INTRODUCTION 

One  of  the  most  outstanding  landscapes  of  grandiose  scale  on  the  earth, 
the  Himalayas  and  the  Qinghai-Xizang  plateau,  has  long  held  the  interest  of 
earth  scientists.  Many  recognized  that  in  this  vast  region,  having  the 
highest  as  well  as  the  youngest  mountain  ranges  in  the  world,  lies  an  impor¬ 
tant  key  to  a  major  tectonic  process  that  took  place  in  the  recent  past. 

Besides  the  surface  expressions,  the  plate  collision  between  India  and 
Eurasia  must  have  produced  fundamental  modifications  in  the  crust  and  upper 
mantle  underlying  this  giant  plateau.  There  is  great  interest  to  delineate 
these  modifications.  A  detailed  delineation  will  undoubtedly  provide  useful 
constraints  on  models  proposed  for  the  collision  process.  Among  geophysical 
tools,  surface  wave  dispersion  analysis  is  one  of  the  most  effective  means 
for  the  crustal  and  upper-mantle  delineations.  Knowledge  of  the  crustal 
thickness  and  velocity  structure  under  the  Qinghai-Xizang  plateau  from  other 
surface  wave  dispersion  studies  is  limited.  Early  works  by  Stonely w  ,  Tandon 
and  Chandhury  C2' ,  Sara  ^  ,  Zeng  et^  aJLkO,  Gupta  and  Narainw,  and  Negi  and 
Singh  O’'  have  given  some  gross  descriptions  of  the  earth's  crust  underneath 
this  plateau.  Most  of  their  results  generally  represent  the  average  crustal 
structure  over  long  surface  wave  paths  through  which  the  plateau  only  accounts 
for  3  j,hort  segment.  More  detailed  works  were  carried  out  by  Bird  and  Toksoz  0/, 
TungLB7f  and  Chun  and  YoshiiC5).  Their  results  were  among  the  first  to  give 
information  on  the  plateau's  layering  crust. 

In  order  to  have  a  surface  wave  study  that  more  directly  reflects  the 
crustal  and  upper  mantle  structure  of  a  tectonic  province,  it  is  important 
to  obtain  as  close  to  pure-path  dispersion  data  as  possible.  Moreover,  re- 
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paated  dispersion  measurements  over  the  same  surface  wave  path  are  desirable 
to  improve  the  data  reliability.  Our  work  begins  by  tectonically  dividing 
China  into  four  major  provinces  here  referred  to  as  subplates  (Fig.  1). 

Each  subplate  has  a  similar  tectonic  past,  therefore,  it  can  be  considered  to 
a  first  approximation  as  a  single  tectonic  unit  over  which  the  lateral  change 
of  crustal  and  upper  mantle  structure  is  relatively  small.  Among  the  four 
subdivisions  in  figure  1,  the  Qinghai-Xizang  plateau  forms  a  subplate  by  it¬ 
self,  and  this  paper  shall  concentrate  on  the  derivation  of  its  crustal  and 
upper  mantle  structure.  To  obtain  pure-path  dispersion  data,  we  take  advan¬ 
tage  of  (1)  several  large  (M  >  7)  earthquakes  that  have  occurred  near  the 
boundaries  of  these  subplates  and  (2)  two  high-quality  Seismological  Research 
Observatory  (SRO)  stations  located  near  boundaries  of  these  subplates. 
Therefore,  a  reasonable  pure-path  configuration  is  closely  approximated.  We 
note  that  tectonically  the  wt-.ern  extension  of  the  Qinghai-Xizang  subplate 
should  also  include  northern  Pakistan  and  Afghanistan,  for  the  entire 
orogeny  belt  was  formed  during  the  collision  process.  Besides  the  two  paths 
crossing  the  Qinghai-Xizang  plateau,  we  elso  studied  three  other  paths  that 
sample  the  surrounding  three  subplates  so  that  a  comparative  analysis  can  be 
made  to  obtain  information  on  lateral  variations  of  crustal  and  upper  mantle 
structure  in  China. 


DATA  AND  METHOD 

Table  1  gives  all  earthquake  sources  used.  They  are  large  aftershocks 
of  the  three  damaging  earthquakes  that  occurred  in  China  in  1976.  A  number  of 
sources  are  used  for  each  path  in  order  to  obtain  an  estimate  on  the  uncer¬ 
tainty  of  the  observed  dispersion  data.  All  events  in  this  study  were  re¬ 
corded  at  two  SRO  stations  located  in  Mashad,  Iran  (MAIO)  and  Taipei,  Taiwan 
(TATO).  SRO  stations  operate  at  high  gain,  thus  lowering  the  detection  thresh 
old  for  seismic  events  and  making  small -magnitude  events  (M  n,  5)  useful  for 
long-  and  short-period  surface  wave  dispersion  analysis.  Using  smaller 
events  would  also  avoid  the  difficulty  in  dispersion  analysis  often  caused  by 
large  multiple  and  propagating  sources.  Since  the  SRO  data  are  digitally 
recorded,  pure  Love  and  Raleigh  waves  are  conveniently  obtained  by  rotating 
the  data  coordinates.  An  example  of  a  seismic  event,  before  and  after 
rotation,  is  given  in  figures  2  and  3.  Surface  wave  data  are  windowed  and 
subjected  to  a  multiple-filtering  analysis.  A  typical  contoured  result  of 
the  relative  amplitude  of  surface  wave  energy  arrivals  over  the  path  Yunnan- 
Mashad  for  both  vertical  and  radial  Rayleigh  waves  is  shown  in  figures  4  and 
5.  The  dots  in  these  figures  represent  the  average  group  velocities  over  the 
filter  band.  Both  fundamental  mode  and  the  first  higher  mode  energies  are 
clearly  present.  Dispersion  data  extracted  from  these  contour  plots  for  the 
two  paths  through  the  Qinghai-Xizang  subplate  are  given  in  figures  6-8. 

Along  the  path  Sichuan-Mashad ,  five  events  each  with  two  (vertical  and  radial) 
Rayleigh  components  are  used  to  obtain  the  observed  dispersion  data  given  in 
figure  9.  Thus  in  figure  9,  every  observed  dispersion  data  point  is  the 
average  of  ten  independent  measurements .  The  standard  errors  as  a  result  of 
these  repeated  measurements  range  from  3-51.  The  averaged  dispersion  data 
for  the  path  Yunnan-Kashad  are  given  in  figure  10.  For  the  purposes  of  com¬ 
parison,  three  other  paths  have  also  been  studied.  Their  averaged  observed 
dispersion  data  are  given  in  figures  11-14.  A  visual  examination  of  these 
observed  dispersion  data  immediately  reveals  that  those  paths  traversing 
through  the  Qinghai-Xizang  plateau  give  substantially  lower  group  velocity 
values.  This  anomalously  low  group  velocity  extends  to  a  period  as  long  as 


80  s,  suggesting  that  the  anomalous  structure  of  the  Qinghai-Xizang  plateau 
extends  down  to  the  upper  mantle. 

A  generalized  linear  inversion  method  first  introduced  by  Backus  and 
Gilbert l1#'is  applied  to  the  observed  dispersion  data.  Detailed  descriptions 
of  this  method  have  been  presented  by  Jacksont11',  WigginsC12\  and  Crosson^13} 
The  basic  logic  of  the  inversion  theory  and  iome  experience  in  its  numerical 
applications  can  be  found  in  Pines  et  a^.Oy.  Initial  models  used  for  the 
inversion  were  based  on  the  results  of  TungC8'.  The  data  were  first  inverted 
for  both  shear  velocity  and  density.  During  the  inversion  process,  large 
instabilities  of  the  density  parameter  were  found  in  the  crust.  This  was 
corrected  by  changing  the  starting  models  to  ones  that  more  closely  approxi¬ 
mated  the  dispersion  data,  and  at  the  same  time  compressional  velocity  was 
also  included  as  an  inversion  parameter^  As  discussed  by  WigginsC12',  Der 
et  aL  Li5),  and  Jackson  and  Burkhardt?  v,  the  partial  derivatives  for  shear 
velocity  are  larger  than  those  for  compressional  velocity  or  density.  There¬ 
fore,  the  resolution  for  shear  velocity  is  considerably  better  than  for  the 
other  two  parameters. 


RESULTS  AND  DISCUSSION 

Upon  application  of  a  generalized  linear  inversion  to  the  observed  dis¬ 
persion  data  (shown  by  open  squares  in  figures  9-14),  we  have  obtained  the 
best-fit  models  corresponding  to  each  path  with  theoretical  group  velocities 
for  each  of  these  models  given  as  open  circles  in  the  respective  figures. 

For  the  two  paths  crossing  the  Qinghai-Xizang  plateau,  the  first  higher  mode 
data  are  available  and  they  were  also  entered  into  the  inversion  process. 

Love  waves  were  not  used  in  the  inversion  due  to  the  absence  of  amplitude 
maxima  in  the  dispersion  curves  between  30  s  and  60  s.  Earlier  studies  by 
Thatcher  and  BruneO7'  and  James C18^  have  discussed  the  problem  of  inter¬ 
ference  of  Love  waves  crossing  inhomogeneous  media.  Tung(0  attributed  this 
problem  on  the  Qinghai-Xizang  plateau  to  Love  waves  following  a  non-least 
time  path.  A  possible  structural  implication  is  that  there  may  be  structures 
along  the  path  corresponding  to  wavelengths  of  Love  waves,  causing  part  of 
the  wave  energy  to  travel  along  a  non-least  time  path  with  slower  velocities. 
Resulting  models  are  labeled  SMI  for  the  Sichuan-Mashad  path,  and  YM1  for  the 
Yunnan-Mashad  path.  Numerical  values  of  the  two  models  are  given  in  Tables 
2  and  3.  In  the  inversion  process,  the  compressional  velocity  and  the  den¬ 
sity  values  were  practically  unchanged  as  they  are  not  very  sensitive  para¬ 
meters  for  the  group  velocity.  Initial  and  final  values  of  the  shear 
velocity  are  tabulated  to  show  the  changes  as  required  by  the  data.  To  com¬ 
pare  the  crustal  and  upper  mantle  structure  of  the  Qinghai-Xizang  plateau 
with  its  surrounding  regions,  we  have  a  summary  plot  given  in  figure  15  where 
the  shear  velocity  structures  of  the  models  SMI  and  YM1  are  shown  together 
with  three  other  models  resulting  from  the  dispersion  data  along  the  three 
other  paths:  TM1  for  Tangshan-Kashad ,  TT1  for  Tangshan-Taipei ,  and  SCI  for 
Sichuan-Taipei .  The  comparison  in  figure  15  shows  a  distinctly  thick  crust 
for  the  Qinghai-Xizang  plateau,  with  lower  crust  extending  to  a  depth  of 
about  70  km.  From  the  width  of  the  resolving  kernels  it  is  only  possible  to 
resolve  layers  of  10-15  km  in  the  upper  crust,  15-20  km  in  the  middle  crust, 
and  25-30  km  in  the  lower  crust.  Resolution  in  the  mantle  is  poor,  a  low- 
velocity  layer  in  the  upper  mantle  can  only  be  inferred.  Additional  data 
for  periods  much  longer  than  120  s  and  good  Love  waves  as  well  as  higher  mode 
data  would  significantly  increase  the  resolving  power  of  the  upper  mantle 
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structure.  The  shear  wave  velocity  structure  for  both  paths  crossing  the 
Tibetan  plateau  is  similar.  A  70  km  thick  crust  may  be  broken  into  four 
layers.  The  first  layer  has  shear  velocities  that  range  between  2.87  and 
3.06  km/s,  suggesting  a  layer  approximately  15  km  thick  of  partly  sedi¬ 
mentary  origin.  Tung(e)  and  Hu  et  ah  discussed  the  existence  of  thick 
sedimentary  layers  in  the  Mount  Oolmo  Langma  region  of  the  Himalayas.  Both 
BirchC20'  and  Simmons C22s reported  shear  velocities  for  sedimentary  rocks  at 
10-  to  15-km  depth  of  approximately  3.00-3.08  km/s.  A  second  layer  has 
shear  wave  velocities  ranging  from  3.20  to  3.40  km/s  and  is  approximately 
15  km  thick,  suggesting  a  layer  of  granitic  or  equivalent  composition.  A 
third  layer  has  shear  velocities  between  3.40  and  3.70  km/s  is  approximately^ 
15-20  km  thick,  and  also  seems  to  consist  of  grarnte.  For  comparison, 

Birch  and  Bancroft^22' ,  Birch  and  Sirmons(21'  reported  shear  velocities 
for  granite  under  3-5  kbars  of  pressure  of  3.45-3.70  km/s.  Simmons  reported 
that  granite  with  a  shear  wave  velocity  of  3.79  km/s  at  pressures  of  10  kbars. 
The  fourth  layer  with  shear  wave  velocities  between  3.85  and  3.95  km/s 
corresponds  to  granite  at  pressures  of  approximately  20  kbars,  which  would 
correspond  to  depths  of  60-70  km. 

The  dispersion  curves  and  shear  wave  velocity  profile  for  the  Tibetan 
plateau  are  in  agreement  with  those  values  presented  by  Tung^  and  Bird  and 
Toksozl7'.  Chun  and  Yoshiii^5',  using  Rayleigh  and  Love  waves  for  many 
different  paths  across  the  Qinghai-Xizang  pleateau,  reported  a  crustal 
thickness  of  approximately  70  km  and  observed  a  velocity  reversal  in  the 
middle  crust.  Owing  to  a  15-  to  20-km  resolution  in  the  middle  crust  it  is 
doubtful  a  velocity  reversal  could  be  observed  for  data  over  the  two  paths 
presented. 

Because  of  the  small  amount  of  information  derived  from  inverting 
Rayleigh  waves  in  the  upper  mantle,  the  above  structural  conclusions  are  only 
tentative.  The  upper  mantle  shear  velocities  are  lower  than  either  the 
Gutenburg  earth  model  or  the  Canadian  Shield.  The  velocity  at  the  top  of  the 
mantle  is  approximately  4.4  km/s,  with  indications  of  a  low-velocity  layer 
at  about  90-km  depth.  Since  the  low-velocity  layer  can  only  be  inferred,  the 
existence  of  extremely  low  velocities  of  4.1  km/s  is  tentative.  Low  shear 
velocities  in  the  upper  mantle  suggests  a  high  degree  of  partial  melting  and 
a  low  resistance  to  deformation,  Anderson  et  alX23^ 
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TABLE  1.  Pertinent  Earthquakes  Information  Used  in  this  Study  (Continued) 


TABLE  2.  Model  Fit  for  Path  Szechwan-Mashad 


TABLE  3.  Model  Fit  for  Path  Yunnan-Mashad 


RELATIVE  AMPLITUDE 


Unrotated  seismograms  for  the  path  Yunnan-Mashad,  July  3,  1976  (from  top  to  bottom 
vertical ,  NS,  and  EW) . 


Rotated  se 1 smoqrams  for  the  path  Yunnan-Mashad ,  July 
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Results  of  multiple  filtering  of  Rayleigh  wave  vertical  component,  Yunnan-Mashad, 
May  31,  1976. 
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•oup  velocity  of  radial  Rayleiqh  component  for  five  events  of  the  path 
chuan-flashad. 
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Group  velocity  of  radial  Rayleigh  component  for  five  events  of  the  path 
Sichuan-Mashad. 
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Rayleigh  wave  group  velocity  for  three  events  of  the  path  Yunnan-Mashad. 


Theoretical  and  observed  Rayleigh  wave  group  velocities  of  the  path  Sichuan-Mashad. 


Theoretical  and  observed  Rayleiqh  wave  group  velocity  of  the  path  Tangshan-Mashad. 
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Theoretical  and  observed  Rayleigh  wave  group  velocities  of  the  path  Tangshan-Tai pei . 
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Theoretical  and  observed  Rayleigh  wave  group  velocities  of  the  path  Sichuan-Taipei . 
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Fig.  14.  Theoretical  and  observed 
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